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Finite volume effects and quenched chiral logarithms 
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We have measured the valence pion mass and the valence chiral condensate on lattice configurations generated 
with and without dynamical fermions. We find that our data and that of others is well represented by a linear 
relationship between ml and the valence quark mass, with a non-zero intercept. For our data, we relate the 
intercept to finite volume effects visible in the valence chiral condensate. We see no evidence for the singular 
behavior expected from quenched chiral logarithms. 
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1. INTRODUCTION 

Over the last few years, analytic calculations 
using chiral perturbation theory have indicated 
pathologies in the quenched approximation jl]^. 
One example is the lack of closed quark loops 
contributing to propagation of the quenched r]', 
which leads to a double pole in the quenched rj' 
propagator. This leads to the following singular 
behavior for the valence pion mass [El: 



ml = (2jum) '+« (A 2 ) 
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For lattice calculations, there is some uncer- 
tainty about the mass m to be used in the above. 
Due to the flavor symmetry breaking of staggered 
fermions at finite lattice spacing, the pions are not 
degenerate in mass in current simulations and the 
quenched 77' may be closer in mass to the heavier 
pions than the lighter ones. Interpreting the m 
above as the bare quark mass assumes the 77' is 
closer in mass to the lightest, Goldstone boson, 
pion §. 

Another consequence of the chiral perturbation 
theory calculations is the absence of quenched chi- 
ral logarithms in decay constants, in particular 
in / w . Because there is a lattice version of the 
Gell-Mann-Oakes-Renner (GMO) relation, this 
means that quenched chiral logarithms in the pion 
mass must be compensated by quenched chiral 
logarithms in the quenched chiral condensate. 
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Recently, there have been reports of simulation 
data consistent with the presence of quenched chi- 
ral logarithms, particularly in To assess 
whether the data clearly reveals quenched chiral 
logarithms, we have investigated three points: 

1. What effects do dynamical fermions have on 
the data? Does the signal persist and is it 
suppressed by the dynamical fermion mass? 

2. Are quenched chiral logarithms evident in 
the valence chiral condensate? 

3. Have all finite volume effects been clearly 
controlled? 

This note will detail some of our results from the 
pursuit of these questions. 

2. RELATING ml and (CO 

Since we will be referring to both quenched and 
dynamical simulations, we will introduce a species 
of staggered fermion, denoted by (, which does 
not appear in the action. This valence fermion 
of mass mQ will be used to construct the valence 
pion propagators and the valence chiral conden- 
sate. The mass of the fermions entering the de- 
terminant will be denoted by m. 

Consider an operator tt(x) which creates a pion. 
We can define a function C(m^) by 



C{m i )=ml(m i ) / d A x{n{x)n(0)) 



(2) 



The angle brackets can represent averaging over 
any ensemble of gauge fields, quenched or dynam- 
ical. Provided the volume is large enough and 
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is much less than the mass of any other state with 
the same quantum numbers, m, in this relation is 
the mass determining the decay of the pion prop- 
agator. 

We now rewrite equation (|J) using the eigen- 
values, A, and the density of eigenvalues per unit 
volume in a background configuration [/, p(X,U), 
of the staggered fermion Dirac operator, D(U). 
We choose our conventions such that D (\ (x) = 
iX(x(x) and J dXp(X,U) = 1. Using the U(l) 
symmetry of staggered fermions, we find that the 
eigenvalues come in pairs, ±A, and 

p(X,fJ,m) 



d 4 x (tt(x)<k(Q)) 



dX 



A 2 



(3) 



where p(A, /3, m) is the ensemble average of 
P (X,U). 

We can also write the valence chiral condensate 
in terms of the eigenvalue density as 

p(X,(3,m) 



dX 



A 2 



We therefore have the relation 

ml(m c ) (CC(™ C )} 



C(m c ) 



(4) 
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If C(m^) is known from either simulations or an- 
alytic arguments, this equation connects m 2 and 
(CC( m c)) including any effects of finite volume or 
quenched chiral logarithms. 

We now focus on equation (^) and make a sim- 
ple model for the effects of finite volume. In 
the continuum, chiral symmetry breaking is as- 
sociated with a non-zero value for p(0), since 
<CC(0)) = 271-/5(0). The major effect of finite vol- 
ume can be postulated to be a minimum eigen- 
value, A m i n , below which the average density of 
eigenvalues vanishes. In particular, consider p(X) 
of the form 

f ° 

where a diverges in the continuum limit due 
to the quadratic divergences in (CC(m^)}. In- 
tegrating equation (H) for our model gives, for 

Amin < TTlQ < Ao, 

Amin , / 

h m^{a 



X < A m in 

Amin < A < A (6) 
A < A 



(CC("i c )} = p° 



(7) 



plus terms of 0(m 2 ) and 0(1/to 2 ). The model 

shows that dependence in (CC( m c)) can ^ e 

easily associated with finite volume effects. 

3. RESULTS FOR m 2 and ((({m c )) 

We will be comparing the results from four 
staggered simulations, two quenched and two un- 
quenched. Figure [l] shows to 2 versus ttiq for these 
four simulations as well as the simulation param- 
eters. The 32 3 x 64 quenched simulation is the 
work of Kim and Sinclair [Q ; the rest were done 
at Columbia. The rise for small is the type 
of behavior expected from equation (|l|). Unfor- 
tunately, this kind of rise has been evident for 
some time in simulations and has been widely at- 
tributed to finite volume effects. In Figure || we 
replot the data, this time without dividing by m^. 
The data is well fit by a straight line, with a non- 
zero intercept. 
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Figure 1. m^/mc, versus rri(. 



The correlated data column of Table |l| gives a 
list of the fit parameters we have obtained for fits 
to m 2 , without the covariance matrix, of the form 

rn 2 = c + cim c (8) 

Of course, one should use the full covariance 
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Figure 2. m 2 . versus ttiq. 



matrix while fitting. We have tried single- 
elimination jackknife fits which simultaneously fit 
m-jr for all values of and an appropriate range 
of times, but we have found such fits to be quite 
unstable given the large dimensionality of the co- 
variance matrix. Instead we have decorrelated 
our data. In particular, we have broken our sim- 
ulations into three or four smaller, equal length 
bins, determined m w for the lightest vtlq from the 
first bin, for the next lightest from the second bin, 
etc. and then done an uncorrelated fit of the form 
given in equation (||) . The results are given in the 
column labeled uncorrelated data in Table |l|. 

We now turn our attention to (££(mf )). Figure 
H is a plot of the valence chiral condensate for the 
four simulations we are considering. Notice the 
log-log scale and the many orders of magnitude 
of ni£ spanned by the Columbia simulations. 

Figure | shows (CCOc)) for the j3 = 5.7, ma = 
0.01 simulation, on a linear plot. The line on the 
graph is a fit of the form 



(CC( m c)) = °o + a i m C _ a -i/ m (. 



(9) 



motivated by our simple model and equation 
(0). For the correlated fit shown, we find a = 
0.0063(3), oi = 2.13(1), <z_i = 8.9(5) x 10~ 7 and 
X 2 /dof = 39(16)/2. The value of x 2 is quite large, 
but the jackknifed error on it is also large. This 



Figure 3. (CC( TO c)) versus m^. The lines on the 
graph just connect the points. The 32 3 x 64 data 
has been scaled to the Columbia normalization. 



generally indicates a poorly resolved small eigen- 
value in the covariance matrix. An uncorrelated 
fit to the data gives a = 0.0061(3), oi = 2.15(2), 
a_i = -8.8(8) x 10~ 7 and x7 dof = 0.8/2. Thus 
we have clear evidence for finite volume effects 
in (CC( m c)) an d ^ey are in agreement with the 
expectations of our simple model. 

4. INTERCEPT FOR ml 

We can now use equation (||) to give us an 
understanding of the origin of the non-zero in- 
tercept, co, found in fitting m\. Assume that 
C(m^) is a non-singular function in the range 
Amin < mc < Xq. Then, using the parameteri- 
zation of equation (J^) gives 



C(0)a. 



C(0) C"(0)a_ 



«o 



m c (10) 



plus terms of 0(m 2 ) and 0(l/m 2 ). This expan- 
sion is only valid in the range where m^ > A m i n , 
which is where most of our data lies. The non- 
zero intercept for m\ is easily tied to the existence 
of a smallest eigenvalue, A m i n in the average eigen- 
value density p. 

One can check that our assumptions about 
C(m^) are valid in this region by computing the 
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Table 1 



Parameters for linear fits of m\ versus tuq. 
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Figure 4. (CC( m c)) versus for the (3 = 5.7, 
ma = 0.01 simulation. 



right-hand side of equation (J5|). We have done 
this and find that C(m^) shows no singular be- 
havior, so the power series expansion of the pro- 
ceeding paragraph appears reasonable. 

Thus we have seen how finite volume effects 
can account for a rise in m^/m^ at (3 = 5.7 and 
ma = 0.01. The other 16 3 x N t runs are con- 
sistent with the same finite volume effect, even 
though they are at much stronger coupling. This 
is reassuring, since the rise in m^/m^ is appar- 
ent for both the quenched and dynamical simu- 
lations. For the 32 3 x 64 simulation, Kim and 
Sinclair have measured the same value for m^ on 
both 24 3 and 32 3 volumes. Thus they interpret 
their rise in m^/m^ as due to a quenched chiral 
log. However, their measurements of (CC( TO c)) 
do not show the logarithmic behavior expected, if 



the rise in m\jmQ is really a quenched chiral log. 
Also, their data is reasonably well represented by 
a fit of the for in equation (||) . 

A possible explanation can come from consid- 
ering the GMO relation for infinite volume. We 
then have 

m\ = C(0) + C'(0)m^ 
m^ ao + airriQ ' 

if C(rn^) is a smooth function of m^. Depend- 
ing on the relative sizes of a±/ao and C"(0)/C(0), 
one could see a rise in m\ jm^ without finite vol- 
ume effects or chiral logarithms, until m^ is very 
small. The relative sizes of ai/ao and C"(0)/C(0) 
are generally [3 dependent, so that even at infinte 
volume m\jmQ could have a different form as a 
function of (3. In addition, if one includes the a_i 
term in equation (|Tl|), there are choices for the 
parameters which can give a decrease in m^/m^ 
for moderate values of m^, and then a rise for 
small m^. Such behavior has been reported by S. 
Gottlieb H at this conference. 

We would like to thank S. Kim and D. Sin- 
clair for a preliminary copy of their paper. We 
also thank Pavlos Vranas and Roy Luo for help- 
ful discussions. 
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